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ABSTRACT 

We present Hubble Space Telescope (HST) Wide Field Gamera 3 (WFC3) observa¬ 
tions of the source and lens stars for planetary microlensing event OGLE-2005-BLG-169, 
which confirm the relative proper motion prediction due to the planetary light curve 
signal observed for this event. This (and the companion Keck result) provide the first 
confirmation of a planetary microlensing signal, for which the deviation was only 2%. 

The follow-up observations determine the flux of the planetary host star in multiple 
passbands and remove light curve model ambiguity caused by sparse sampling of part 
of the light curve. This leads to a precise determination of the properties of the OGLE- 
2005-BLG-169Lb planetary system. Combining the constraints from the microlensing 
light curve with the photometry and astrometry of the HST/WFC3 data, we find star 
and planet masses of M* = 0.69 ± 0.02Mq and nip = 14.1 ± O.9M0. The planetary mi¬ 
crolens system is located toward the Galactic bulge at a distance of Dl = 4.1 ± 0.4 kpc 
and the projected star-planet separation is a_L = 3.5 ±0.3 AU, corresponding to a semi¬ 
major axis of a = 4.0(’jpg AU. 

Subject headings: gravitational lensing: micro, planetary systems 


1. Introduction 


Gravitational microlensing is unique among planet detection methods (Bennett 2008; Gaudi 


2012) in its sensitivity to planets with masses smaller than Earth (Bennett & Rhie 1996) orbiting 
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beyond the snow line (Mao k, Paczyhski 1991; Gould & Loeb 1992), where planet formation is 


thought to be the most efficient (Ida & Lin 2005 Lecar 2006; Kennedy et al. 2006; Kennedy k 


Kenyon 

2008 

Thommes et al. 

2008 

), according to the core accretion theory of planet formation 

(]Lissauer 

1993 

Pollack et al. 

1996 

)• 

Microlensing is also able to detect planets orbiting stars at 


distances ranging from a few hundred parsecs up to ~ 8kpc. Since the microlensing method 
doesn’t depend on light from the planetary host star, it can be used to find planets orbiting very 


faint star or even stellar remnants or brown dwarfs (Bennett 2008; Gaudi 2012). However, one 


drawback of the microlensing method is that the microlensing light curves usually do not indicate 
the planet or host star mass. Instead, they generally yield the planet-star mass ratio, g, and the 
separation in units of the Einstein radius (Re), except for events that exhibit the microlensing 
parallax effect ( Gaudi et al.||2008 Bennett et al.|[2008 Muraki et al. 2011 Furusawa et al.|[2013 ). A 
measurement of the microlensing parallax effect for a planetary microlensing event usually provides 
enough information about the lensing geometry to determine the lens mass. The mass measurement 
does require that the angular Einstein radius, Be = Re/Dl, be known, but this can be determined 
for most planetary events from finite source effects in the light curve that allow the source radius 
crossing time, t*, to be measured. However, most events do not have a measurable microlensing 
parallax effect, particularly those due to lens systems in the Galactic bulge. 

A more generally applicable method to determine the lens system mass is to detect the host 
star, for an event in which Be has been determined. This requires high angular resolution imaging 
because the lens and source stars are not resolved from unrelated stars in ground-based, seeing- 
limited images. When Be is known, it provides a mass-distance relation for the lens system, and 


this can be combined with a mass-luminosity relation to determine the mass 

of the lens system. 

This has been done for a number of events (Bennett et al. 

2006 

Dong et al. 

2009 

Kubas et al. 


2012 

Batista et al. 

2014 

), but sometimes it isn’t clear if the excess flux is really due to the lens 

star ( 

Sumi et al. 

2010 

Janczak et al. 2010 

Gould 

2014 

), as unrelated stars or companions to the 


source or lens star cannot always be excluded. The keys to establishing that the excess flux is due 
to the planetary host (and lens) star are to measure lens brightness in multiple pass bands and to 
measure the relative lens-source proper, which is usually known from the light curve. 

In this paper, we present the first direct measurements of the relative proper motion, 
for a planetary microlensing event, OGLE-2005-BLG-169, using HST observations in three Wide 
Field Camera 3 (WFC3) passbands: F814W, F555W, and F438W. The light curve prediction of 
//rel comes from the planetary signal itself, so our confirmation of this prediction is a confirmation 
of the planetary signal. Thus, the planetary signal for OGLE-2005-BLG-169Lb is the first to 
be confirmed by follow-up observations. The HST follow-up observations also provide a tighter 
constraint on /Xrei than the light curve does, so we are able to obtain tighter constraints on the 


light curve parameters than the discovery paper (Gould et al. 2006). The HST lens brightness 


measurements, when combined with the Be mass-distance relation, yield the masses and distance 


of the planet and its host star, as well as their projected separation. A companion paper (Batista 


et al.|20l5) presents independent measurements of /irei the lens brightness in the Lf-band using 






















































































adaptive optics observations from the Keck-Il Telescope. These Keck measurements are consistent 
with the HST results presented here. 


This paper is organized as follows. We discuss the light curve data and photometry in Sectionj^ 
and in Sectionj^we present the light curve models that are consistent with the data. In Section]^ we 
show how the angular radius of the source star relates to its color and brightness. Then in Sectionj^ 
we describe the HST data and its reduction, and in Section 5.1 we compare the lens-source relative 
proper motion prediction from the light curve with the HST measurement. In Section 5.2 


we 


compare our results to the Keck adaptive optics observations made 1.74 years later and show that 
the combined HST and Keck observations conhrm that our identification of the lens star is correct. 
The constraints on the lens system from the HST data are explored in Section 5.3 Finally in 


Section we present our conclusions and explain how this analysis demonstrates the primary 


exoplanet host mass measurement method for the WEIRST and EUCLID missions ( 

Bennett & 

Rhie||2002 

Bennett et al. 

2007; 

Green et al. 

2012 

Spergel et al.||2013 

Penny et al.||2013 

)• 


2. Light Curve Data and Photometry 


OGLE-2005-BLG-169 is unique among planetary microlensing events in a number of respects 
( ]Gould et al. 2006). It has the smallest impact parameter, uq, of any planetary microlensing 
event, and it has the smallest amplitude photometric signal of any planetary microlensing event. 
The planetary signal entirely in the extremely high cadence data taken from the 2.4m MDM 
telescope. (More than 1000 observations were taken in a 3-hour period at high magnification.) 
Because of the low amplitude signal, there was concern that the data could be contaminated 
by systematic photometry errors. Due to this concern, the MDM data were reduced with two 


independent photometry pipelines, the OGLE pipeline (Udalski 2003) and the Hartman et al 


(2004) implementation of the Alard & Lupton (1998) photometry code. This later reduction was 


performed by K.Z. Stanek, and we will refer to it as the Stanek reduction. 

In addition to the MDM data set, which contains the planetary signal, the photometry for this 
event include data from the 1.3 m OGLE survey telescope (responsible for the identification of the 
microlensing event), the 1.3m SMARTS telescope at GTIO in Chile, the 2.0m Eaulkes Telescope 
North in Hawaii, and the 0.35 m Nustrini Telescope in Auckland, New Zealand. We use the same 


photometric reduction for each data set that was used by Gould et al. (2006) except for the GTIO 


data. A minor problem was discovered in the GTIO /-band photometry used in the original paper. 
The GTIO /-band photometry yielded a source magnitude that was 0.13 mag fainter than the 
source magnitude from the OGLE /-band photometry when the photometry from both data sets 


was calibrated to the OGLE-HI photometry database (Szymahski et al. 

2011 

). This inconsistency 

was largely resolved (reduced to 0.03 mag) by switching from DoPHOT ( 

Schechter, Mateo, &: Saha 


GTIO //-band data, taken simultaneously with the V and /-band data on the Andicam instrument 
on the SMARTS telescope. The //-band data is especially useful because they allow a more precise 
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3491.8 3491.9 3492 

HJD - 2450000 


Fig. 1.— The light curve peak of event OGLE-2005-BLG-169 with photometric measurements from 
the MDM 2.4m /-band (red, Stanek reduction), OGLE /-band (black), CTIO I and //-bands (blue 
and green), Eaulkes Telescope North (magenta), and the Auckland unfiltered telescope (gold). The 
best ht model is indicated by the black curve, and the grey dashed curve indicates the same model 
without the planetary signal. The bottom panel shows the residual with respect to this no-planet 
model. The MDM data clearly trace out the caustic exit feature, but the data on the rising side 
provide a very weak constraint on the caustic entry properties. So, a variety of angles between 
the lens axis and the source trajectory are permitted by the photometry. The best fit model to 
the data set including the Stanek MDM photometry presented here is consistent with our proper 
motion measurement, while the light curve presented in Gould et al. (2006) is not. 
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determination of the angular source radius (Kervella et al. 2004; Boyajian et al. 2014). The H- 
band light curve used in this paper is a SoDoPHOT reduction, but a reduction using the MOA 
Collaboration difference imaging pipeline (Bond et al. 2001) gives indistinguishable results. 


3. Light Curve Models 


The light curve models used for this paper are different from the models presented in Gould et 
al. (2006) because a different data set is used. We use the Bennett|(2010) modeling code instead of 


the Gould et al. (2006) code, but this has no effect on the results, as these codes have been shown to 
give identical results to better than 1 part in 10^. Our conclusions based on the light curve modeling 


alone are essentially the same as the conclusions of Gould et al. (|2006) . As discussed in Gould 


et al. (2006) and Section]^ the planetary signal for this event is particularly sensitive to potential 
systematic photometry errors because of the larger than usual S/N of the MDM observations and 


the small amplitude of the planetary signal. For this reason, Gould et al. (2006) did the complete 


analysis using both the Stanek and OGLE-pipeline photometry. We continue this philosophy in this 
paper and assume that the OGLE-pipeline and Stanek reductions are equally likely to be correct, 
and so we perform Markov Ghain Monte Garlo (MCMG) calculations for both data sets starting at 


the parameters of each of the local minima presented in Gould et al. (2006). 


The results of these MGMG calculations differ in detail from the results presented in [Gould et 


al. (2006), in the sense that the models with the source trajectory nearly perpendicular to the lens 


axis are now somewhat favored with respect to the previous analysis. With the Stanek version of 
the MDM photometry, these models are now favored by Ax‘^ = 8.8 over the best fit model with a 
source trajectory > 25° from perpendicular to the lens axis. The best fit model using the Stanek 
version of the MDM photometry is presented in Figure and the parameters of this model and 
the best fit s < 1 model are given in Table The best fit model, which has s > 1, is labeled as 
“Stanek s > 1” and the parameters of the best fit s < 1 are also given. This s < 1 model is very 
slightly disfavored with = 0.12. Table also gives the MGMG averages of the parameters 
both without the constraints from the HST measurements (in the next-to-last column) and with 
the constraints from the HST measurements in the last column. Because of the wide variation in 
the 6 values (source trajectory angles) allowed by the light curve, there is a large scatter in some 
of the other fit parameters, such as the source radius crossing time, t*, and the planet:star mass 
ratio, q. 

Figure shows a close-up of the caustic configuration for the best-fit model with the source 
trajectory given by the solid black line. The red circle indicates the size of the source star, and the 


gray dashed line shows the source trajectory for the model presented in Gould et al. (2006). 
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Fig. 2.— The caustic configuration for the OGLE-2005-BLG-169 model shown in Figure The 
black line, with arrow, shows the source trajectory for this model, while the grey dashed line shows 
the source trajectory for the other local minimum for the light curve modeling. This and similar 
models are consistent with the light curve, but they are contradicted by the relative proper motion 
measurement that we present here. 
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Fig. 3.— The {V — /, I) color magnitude diagram (CMD) of the stars in the OGLE-III catalog 
( jSzymahski et ah 2011) within 120" of OGLE-2005-BLG-169. The red spot indicates red clump 
giant centroid, and the blue spot indicates the source magnitude and color. 








4. Calibration and Source Radius 


In order to measure the angular Einstein radius, Oe = we must determine the angular 

radius of the source star, 6^, from the dereddened brightness and color of the source star ( jKervella 
et al.|2004 Boyajian et al.|2014 ). We determine the source star brightness in the V and /-bands by 


calibrating the CTIO E-band and OGLE /-band magnitudes to the OGLE-III catalog (Szymahski 


et al. 2011) yielding the following relations: 


Vbscai = 23.08516 0.97257 Esod + 0.02743 /qsic ± 0.004 

-^oscai = 1.406255 0.93933 /o3ic + 0.060674 Esod ± 0.004 . 


( 1 ) 

( 2 ) 


/o 3 ic is the OGLE /-band light curve magnitude, which differs from the standard Cousins /-band 
used in the OGLE-III catalog, while Eoscai and /oscai refer to the Johnson E-band and Cousins 
/-band magnitudes, as presented in the OGLE-III catalog ( Szymahski et al.]|2011 ). Esod is the raw 
CTIO E-band photometry from our SoDoPHOT reduction. The E-band calibration is based on 54 
stars with 1.0 < (E — /)o 3 cai < 2.2 and /o 3 cai < 16.0 within 2 arc minutes of the target star, and 


the /-band calibration employs the formulae presented in Szymahski et al. (2011). 


Our CTIO //-band SoDoPHOT magnitudes are calibrated to 2MASS ( |Carpenter||2001[ ) with 
the following relation, 

//2mass = //sod + 19.849 ± 0.010 , (3) 

based on 36 stars within 105'' of the target. 


Table 1. Model Parameters 


parameter 

units 

Stanek s > 1 

Stanek s < 1 

MCMC 

no const. 

averages 

//pel H const. 

tE 

days 

43.09 

43.16 

41.8(2.9) 

42.5(1.4) 

to 

HJD - 2453490 

1.8784 

1.8784 

1.8776(10) 

1.8784(1) 

^min 


0.001229 

0.001228 

0.001267(9) 

0.001250(4) 

s 


1.0190 

0.9828 

1.004(18) 

1.001(18) 

9 

radians 

1.6025 

1.6069 

1.43(20) 

1.60(3) 

q 

10-5 

5.913 

5.844 

7.07(1.22) 

6.15(30) 


days 

0.02174 

0.02168 

0.0202(17) 

0.0228(5) 

Oe 

mas 

0.905 

0.911 

0.965(94) 

0.848(27) 

ltrel,G 

mas/yr 

7.67 

7.69 

8.47(87) 

7.29(15) 

Hs 


18.852 

18.854 

18.81(8) 

18.84(4) 

Is 


20.592 

20.594 

20.55(8) 

20.58(4) 

Vs 


22.254 

22.257 

22.21(8) 

22.24(4) 

fit 


1146.66 

1146.78 
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In order to estimate the source radius, we need extinction-corrected magnitudes, and we de¬ 
termine these from the magnitudes and colors of the centroid of the red clump giant feature in the 
color magnitude diagram (CMD), as indicated in Figure]^ The extinction can be determined most 
accurately if three colors are used ( Bennett et al.|[2010[ ), and we find that the red clump centroid 
in this held is at Id = 15.61, {V — I)c\ = 1.93, (/ — H)c\ = 2.07, which implies Hc\ = 13.54 and 
Ki = 17.54. 


We follow the method of Bennett et al. (2010) to determine the extinction, but we use the 


updated dereddened red clump magnitudes of Nataf et al. (2013). We assume absolute red clump 
giant centroid magnitudes of = —1.30, M/d = —0.13, and My^x = 0.93. The Galactic 

coordinates of OGLE-2005-BLG-169 are ((,5) = (0.6769°,—4.7402°), and this implies a distance 
modulus of DM = 14.541. Using the Bennett et al. (2010) method, we estimate the extinction 


toward the center of the Galaxy in this direction to be Ah = 0.374 ± 0.020, Aj = 1.256 ± 0.050, 
and Ay = 2.132 ± 0.090. These extinction values allow us to determine the dereddened magnitude 
for each passband, Cso = Cg — Ac where C refers to the passband (either V, I, or H). 


These dereddened magnitudes can be used to determine the angular source radius, 0*. Of the 
measured source magnitudes, the most precise determination of 0* comes from the (V — H),H 
relation. We use 


log/o [20*/(lmas)] = 0.536654 -h 0.072703 (V - H)go - 0.2 Hgo , 


(4) 


which comes from the Boyajian et al. (2014) analysis. These numbers are not included in the 


Boyajian et al. (2014) paper, but they were provided in a private communication from T.S. Boyajian 


(2014). She reports that this formula determines 0* better than 2% accuracy. This is somewhat 
better than the 2.6% accuracy of the {V — H), H relation of Kervella et al.| ( 2004[ ). (They report 
an accuracy of 1.12% for log/o(0*), which corresponds to 2.6% accuracy for 0*.) 


The implied source radii for the best fit s > 1 and s < 1 models are given in Table along 
with the angular Einstein radius, 0/;, and the lens-source relative proper motion, Hrei,G = 0*1^*, in 
a geocentric reference frame. The light curve parameter that the relative proper motion depends on 
is the source radius crossing time, t*, which is measured in the reference frame of the Earth-bound 
observatories that observe the light curve. So, t* is measured the Geocentric reference frame moving 
at the instantaneous velocity of the Earth at the time of the event, and this is the reference frame 
that /irei,G is determined in. The green histogram in Eigure|^ shows the distribution of Hrei,G from 
our MCMG light curve modeling calculations using both the Stanek and OGLE-pipeline reductions 
of the MDM data. The spread in /irei,G values is primarily due to the uncertainty in the source 
trajectory angle, 0, as discussed in Section]^ 

In Section we will present the relative lens-source proper motion measurement from the 
HST observations. This measurement is made with respect to the average motion of the Earth 
during the 6.4678 years between the event and the HST observations. If we assume that the HST 
observations are made in a Heliocentric frame, the maximum error in the lens-source displacement 
is twice the relative lens-source relative parallax or 27ri.ei = 2A\J{1/Dl — 1/Ds) — 0.26 (assuming 
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6 8 10 12 
Mrei (Geocentric) 

Fig. 4.— Comparison of geocentric proper motion predictions from all Galactic bulge microlensing 
events (grey histogram), microlensing models consistent with the light curve data (green histogram) 
and with the HST measurement (red cross-hatched histogram). The light curve /irei,G distribution 
is drawn from MCMC calculations using both the Stanek and OGLE MDM reductions. The 
HST /irei,G distribution has been converted from the two-dimensional / 2 j.g[ measurement using a 
probability distribution for relative lens distance (Dl/Ds) from a Galactic model. 
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our final result), which compares to our lens-source displacement measurement error of 1.3mas, so 
the assumption of a Heliocentric reference frame is a reasonable approximation. (The lens-source 
relative parallax is given by TTrei = AU ) The HST measurements also determine the 

direction of the lens-source relative proper motion, so they determine the 2-dimensional relative 
proper motion, /Xj.gj jj. 


5. HST Astrometry and Photometry 


We observed the OGLE-2005-BLG-169 source and lens stars for two HST orbits as a part of 
HST Program GO-12541. On 2011 October 19, we obtained images in three passbands, F814W, 
F555W, and F438W, using the Wide Field Camera 3-Ultraviolet-Visible (WFC3-UV) instrument. 
We obtained 7 x 85 sec dithered F814W exposures, 8 x 175 sec dithered F555W exposures, and 
6 X 349 sec dithered F438W exposures. Due to the relatively short exposure times, we have been 
forced to limit the amount of data read out for the F555W and F814W exposures. Only Hex Ik 
region of the CCDs were read out for these passbands. 


The data were reduced following the method of Anderson & King (2000, 2004). The dithered 
exposures are used to construct an effective PSF from stars of a similar color to the target (i.e. 
the blended image of the source plus lens stars). Then, this effective PSF is used to fit two stellar 
profiles to the blended target image. The top-right and bottom two panels of Figure show close- 
ups of the blended source plus lens stars in the three passbands, F814W, F555W, and F438W, 
which are the HST versions of the I, V, and H-bands. The best fit locations of the lens and source 
stars are also indicated. In the F814W images, both the stars have a brightness consistent with 
/-band source brightness determined from light curve modeling, so there would be ambiguity in 
the lens and source star identifications if we had data in this passband alone. Fortunately, the lens 
is considerably fainter than the source in the F555W and F438W passbands, and this allows us to 
uniquely identify the lens and source stars. (The lens is closer than the main sequence source, so 
it must be redder than the source if it has the same magnitude in the /-band.) 


The separation between the lens and source stars is due to the At = 6.4678 yr interval between 
the event and the HST observations, so we can determine the relative proper motion in the Helio¬ 
centric frame by /.treiH = Ax/At, where Ax is the two-dimensional separation between the lens 
and source stars as measured in the HST images. In the Galactic coordinate system, we find 


= (7.52 ±0.27,1.07 ±0.28) mas/yr (F814W) , (5) 

= (7.17 ± 0.33,1.88 ± 0.41) mas/yr (F555W) , (6) 

Hrei,H(^^) = (7.32 ±0.67,1.40 ±0.82) mas/yr (F438W) . (7) 


The error bars for these / 2 j.gi values are determined from the dual star fits. First the fit values 
are renormalized to give y^/d.o.f. = 1. The original y^/d.o.f. values were 1.25, 1.64, and 1.39 for the 
F814W, F555W, and F438W passbands, respectively. Then, we add 0.1 mas/yr in quadrature and 
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Fig. 5.— The top-left panel shows a 4.9" x 4.6" section of the stacked F814W images containing 
the target, and the top-right panel shows a close-up of a 0.42" x 0.37" region containing the lens 
and source stars. The lower-left and lower-right panels show the same region using the stacked 
images in the F555W and F438W passbands. The magenta spots on the left are the best fit lens 
locations and the cyan spots on the right are the best fit source positions. These close-up images 
are sums of all the dithered images with 100 x oversampling. The distortion of the WFC3 images 
has been removed, and this results in pixels shaped like parallelograms rather than squares. 
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multiply the error bars by 1.5. These adjustments are meant to account for systematic uncertainties 
in the PSF models, and they ensure that values from the different passbands are consistent. (The 
systematic errors could probably be reduced by constructing different PSF models for the lens and 
source stars rather than one PSF based on their average color.) Note that the F438W and the Keck 


Ff-band (Batista et al. 2015) proper motion values both fall in between the F814W and F555W 


values, so there is no trend with color. We combine the best two measurements from equations [5]j^ 
(F814W and F555W), to obtain the measurement we will use as our final measurement of the 
lens-source relative proper motion 


/^rei,H(^ = (7.39 ± 0.20,1.33 ± 0.23) mas/yr . 


( 8 ) 


The direction of proper motion is about ~ 10° from the Galactic longitude, I, direction. This 
is expected for a Galactic disk lens about half way to the center of the Galaxy. Due to the Galaxy’s 
flat rotation curve, we and the lens system move at about the same velocity, but the source star in 
the Galactic bulge doesn’t share this rotation, so the relative lens-source proper motion is typically 
Hrei ~ 220 km/s/8.3 kpc = 5.6 mas in the direction of Galactic rotation. The velocity dispersion is 
dominated by the Galactic bulge one-dimensional velocity dispersion of ~ lOOkm/s, which implies 
a one-dimensional proper motion dispersion of 2.5 mas/yr. So typically, the relative proper motion 
of lens half way to the bulge should be within ~ 30° of the Galactic rotation direction, while an 
event like OGLE-2005-BLG-169, with a higher than average /Xreb would typically have a relative 
proper motion within ~ 30° of the Galactic rotation direction. 


These fits also return the magnitudes of the source and lens stars. To put these on a standard 
scale, we calibrate the V (F555W) and I (F814W) to the OGLE-III catalog ( Szymahski et al.| 2011). 
We find 12 uncrowded stars with 0.86 < V — I < 2.24 and / < 19 that we use for this calibration. 
The scatter in these calibrations are about 1%, so the formal uncertainty in the calibration is 
0.3-0.4%. 


With this calibration, the best fit V and I source magnitudes are 

(9) 
( 10 ) 
( 11 ) 


Ks = 22.212 ± 0.041 Is = 20.555 ± 0.054 

Vl = 22.783 ± 0.067 II = 20.493 ± 0.051 

= 21.704 ± 0.020 Itot = 19.771 ± 0.020, 


where Ftot and /tot refer to the magnitudes corresponding to the combined brightness of the lens 
plus source stars. The formal errors on 14ot and Itot are actually only about 0.003 mag, but we 
use 0.020 mag to account for calibration uncertainties. The uncertainties on the source and lens 
magnitudes are significantly larger than the uncertainty on the combined lens-|-source magnitude. 
This is due to the fact that the lens and source are not fully resolved, which allows correlated 
uncertainties where the lens and source can trade flux with slight modifications in their best fit 
positions. 

The E438W data were calibrated to the Vega-magnitude scale by comparison to a reduction 
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of the same data using DolPHOT, which is an updated version of HSTphot ( Dolphiri||2000 ) by the 
same author. This gives Bs = 23.82 ± 0.06 and = 24.74 ± 0.15. 


5.1. Confirmation of Planetary Signal Prediction 


Our HST measurements of the lens-source relative proper motion, /Xj-guj, are made in a ref¬ 
erence frame that is indistinguishable from the Heliocentric reference frame, but the light curve 
measurements are made in the Geocentric reference frame that moves with the Earth at the time 
of the event. These reference frames differ by the velocity of the Earth at the time of the event 
projected onto the plane of the sky at the time of the event. This projected velocity is 

{v!^n,v®e) = (3.15,18.51) km/sec = (0.665,3.905) AU/yr , (12) 


and the relationship between the geocentric and heliocentric relative proper motions is 


TTf-gl 

Prel.H = A'rel.G + 


(13) 


Converting to Galactic coordinates, we have (v^ijV^b) = (3.74,—1.30) AU/yr, so equation 13 be¬ 
comes 

(14) 


Erel , 




in Galactic coordinates or 


Prel,G(^^) — (7.39,1.33) -|- 


■^rel 


(-0.49,0.17) , 


(15) 


0.13 mas 

after substituting our measured value from equation!^ We choose vTrei = 0.13mas as our reference 


value in equation 15 because this is a round number close to our best fit final value. It corresponds 
to a lens system at a distance of about 4kpc, about half-way to the Galactic center. 

However, we’d like to compare our measurement of /Xj-gj jj to the light curve prediction of jiirei,G 
with no reference to our best fit lens distance, in order to have a relatively pure test of the relative 
proper motion prediction from the light curve. We therefore convert the measured Heliocentric 
relative proper motion measurement to a probability distribution for /Xrei.G using a Bayesian analysis 


with the Galactic model of Bennett et al. (2014). This analysis implicitly assumes that potential 


primary lens mass for the OGLE-2005-BLG-169 event is equally likely to host a planet with the 
measured mass ratio, but it makes no assumptions about the location of the lens system or source 
star. 


The red cross-hatched histogram in Figure [^indicates the distribution of ^rei,G values consistent 
with the HST /Xj.gijj measurement. The HST values for jiirei,G are near the extreme low-/rrei,G edge 
of the light curve distribution, shown in green. However, the histograms cross at a probability of 
about 80% of the maximum of each curve. Thus, the HST measurement is clearly consistent with 
the /irei,G predictions from the light curve. 
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The light curve prediction of /Trei.G comes directly from the planetary light curve feature, 
because this is the only light curve feature that resolves the finite source size and determines the 
source radius crossing time, t*. Thus, our HST confirmation of the /irei,G is also a confirmation 
of the planetary interpretation of the OGLE-2005-BLG-169 light curve. This is the first such 
confirmation for a planetary microlensing event. 


5.2. Comparison to Keck Adaptive Optics Measurements 


In July, 2013, a subset of us obtained 15 Keck NIRC2 Adaptive optics ff-band images with 
seeing of ~ 55 mas. These high resolution images resolved the lens and source stars, so that their 
separation could be measured, some 8.2121 years after the microlensing event peak. This allowed 
an independent measurement of the lens-source relative proper motion (Batista et al. 2015), 


Pj.ei h(^ &)[Keck] = (7.28 ± 0.12,1.54 ± 0.12) mas/yr . (16) 

This measurement is obviously quite consistent with our HST measurement given above (equation]^ 
as both the I and b components are within l-a of our values. 

Both of these ^ measurements make the assumption that the lens and source are coincident 
during the microlensing event, but we can also use the 1.7443 yr interval between the HST and Keck 
observations to work out the separation at the time of the event between the stars we identify as 
the lens and source. This gives a separation of 


(Al, Ab) = (3.5 ± 7.2, -6.4 ± 7.8) mas , 


(17) 


between these lens and source stars at the time of the event. These are consistent at < Icr with 
our identification of the lens and source stars, whose separation was 9F,U m\n = 0.0011 mas at the 
time of the event. 


This measurement is also sufficient to rule out the possibility that the detected flux is due 
to a binary companion to the lens. A possible binary companion to the lens that orbits within 
30 mas of the primary lens star would strongly perturb the light curve, so such a close companion 


is excluded by the light curve observations, which see no evidence of such a companion (Bennett et 


al. 2007 Gould 2014). Thus, this light curve constraint combined with the combined Keck+HST 


measurements, excludes the possibility that a binary companion to the lens is responsible for the 
flux that we attribute to the lens stars. 


5.3. Lens System Properties from HST Measurements 

As discussed in Section]^ the angular Einstein radius, 6e = O^tE/t^, can be determined from 
light curve parameters, as long as the angular source size, 0*, can be determined from the source 
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brightness and color. The determination oi 9 e allows us to use the following relation (Bennett 


2008; Gaudi 2012 ) 




9e 

1 mas 


I — X 


Ds 

S kpc 


(18) 


where x = Dl/Ds- This expression can be considered to be a mass-distance relation, since Ds is 
approximately known. As can be seen from Table the light curve does not determine 0e very 
precisely, due to the correlated uncertainty in 9 and t*. However, as Figure indicates, the HST 
observations rule out a large fraction of the /irei.G values that are compatible with the light curve. 
Since 9e = this implies that much of the 9e range allowed by the light curve is now 

excluded by the HST data. This, in turn, has an effect on other parameters, such as the planet:star 
mass ratio, which is q ~ 6 x 10“^ for the 9 ~ 1.6 solutions, compared to (7 ~ 8 x 10“® for the 9 ~ 1.0 
solutions. So, the HST data drive the planetary mass fraction to a somewhat lower value. 


To solve for the planetary system parameters, we sum over our MCMC results as in Subsec¬ 
tion 5.1 with the weighting by the Galactic model parameters consistent with the HST /x^ein mea¬ 
surement, but now we add the HST lens brightness constraints, as well. In this sum, we randomly 
select source and lens distances that are consistent with the mass-distance relation (equation |18[). 
In order to check this consistency, we must invoke a mass-luminosity relation. We use the mass- 
luminsity relations of Henry & McGarthy (1993), Henry et al. (1999) and Delfosse et al. (2000). 
For Ml > O.Q6Mq, we use the Henry & McCarthy (1993) relation; for 0.12 Mq < Ml < 0.54M©, 


we use the Delfosse et al. (2000) relation; and for 0.07 Mq < Ml < 0.10 Mq, we use the Henry et 


al. (1999) relation. In between these mass ranges, we linearly interpolate between the two relations 


used on the boundaries. That is we interpolate between the Henry & McCarthy (1993) and the 


Delfosse et al. 

to 

0 

0 

0 

) relations for 0.54 Mq 

< Ml < 0.66Mq, and we interpolate between the 

Delfosse et al. ( 

to 

0 

0 

0 

and 

Henry et al. 

dlMo 

) relations for 0.10 Mq < Ml < 0.12 Mq. 


Table 2. Physical Parameters 


Parameter 

units 

value 

2-a range 

Dl 

kpc 

4.1 ±0.4 

3.3-4.8 

M* 

Mq 

0.69 ±0.02 

0.64-0.73 

rrip 

Mq 

14.1 ±0.9 

12.4-15.9 

a± 

AU 

3.5 ±0.3 

2.9-4.0 

asd 

AU 

4 

^•'^- 0.6 

3.0-14.0 


Note. 

ranges. 


Uncertainties are 1-a parameter 
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At a Galactic latitude of 6 = —4.7402°, and a lens distance of ~ 4kpc, the lens system is likely 
to be behind most, but not all, of the dust that is in the foreground of the source. We assume a 
dust scale height of /ijust = 0.10 ± 0.02kpc, so that the extinction in the foreground of the lens is 
given by 


I _ g-|.Di/(/iduat sinb)! 

1 — e-|Os/(/idustsinfe)| > 


(19) 


where the index i refers to the passband; V, I, or H. For each model in the Markov Chain, the 
^dust value is selected randomly from a Gaussian distribution. We assume error bars of ay = 0.10 
and aj = 0.07 magnitudes for the combined uncertainty in the mass-luminosity relations and 
the lens star extinction estimate. The results of this hnal sum over the Markov Chain are given 
in Table The host star is a M* = 0.69 ± 0.02 Mq K-dwarf, orbited by a planet of about 
Uranus’ mass at rup = 14.1 ± 0.9M 0 , at a projected separation of a± = 3.5 ± 0.3 AU. Assuming 
a random orientation, this implies 3-dimensional separation of a^d = 4.0 _gg AU. This planet then 
has the mass of an ice-giant in a Jupiter-like orbit at about twice the nominal snow-line distance 
of 2.7{M/Mq) AU ~ 1.9 AU. This is similar to a number of other planets found by microlensing 
( ]Beaulieu et al. 2006} Sumi et al. 2010 Muraki et al. 2011, Furusawa et al. 2013), which can be 
interpreted as examples of “failed Jupiter cores’.’ These would be planets that grew by accumulation 
of solids, as Jupiter’s core is thought to have done ( jLissauer 1993), if the core accretion model is 
correct. These “failed Jupiter core” planets are thought to be common around the low-mass stars 
probed by the microlensing exoplanet search method. 


6. Discussion and Conclusions 


We report the first detection of an exoplanet microlens host star found at the separation 
predicted by the exoplanet feature in the microlensing light curve. Together with a companion 
paper based on Keck data ( Batista et al.|20T5 ), this provides the first confirmation of a microlensing 
planetary signal, in the sense that the planetary interpretation of a light curve feature predicted 


the lens-source relative proper motion, which we, and Batista et al. (2015) have confirmed. 


The resulting system is a Uranus-mass planet orbiting a K-dwarf at about twice the snow-line. 


which hts the properties of a “failed Jupiter core” planet, predicted by core accretion (Laughlin et 


al. 2004). 


This is also the hrst demonstration of the primary exoplanet host star mass measurement 
method (Bennett et al. 2007) planned for WFIRST ( jSpergel et al. 2013) and EUCLID (Penny et 
al. 2013). While, the host star mass might plausibly be inferred from just the brightness of the 


lens star (Bennett et al. 2006; Dong et al. 2009 Batista et al. 2014), but to be highly confident 


that the measured star is actually the lens star (Janczak et al. pMol ), it is necessary to measure the 
lens-source relative proper motion and show that it is consistent with the prediction from the light 
curve. This will be even more important for the WFIRST exoplanet microlensing survey, because 
it will work in more crowded helds, where the microlensing rate is highest. 
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The HST data presented her provide an extremely high S/N measurement of the lens-source 
relative proper motion. The lens-source separation is measured at 28-a in the F814W band, 22-a 
in the F555W band, and 11-a in the F438W band. This is partly because this event is a favorable 
one for such measurement (Henderson et al. 2014), but also because it took a while for the HST 
TAG to recognize the importance of such measurements. As discussed in Bennett et al. (2007), this 
measurement could easily have been made 4 years earlier. 


It was not necessary to achieve the photon noise limit in our HST astrometry measurements 
because of the high S/N in the HST data. As the discussion in Section [^indicates, our error bars 
are probably about a factor of two above the photon noise limit. There are several things that 
can be done to improve the analysis. One improvement would be to add a second iteration of PSF 
fitting to determine the source and lens properties. The first iteration determines the approximate 
source and lens colors, but the stars selected to make the PSF models are matched to the average 
lens-l-source color. In a second iteration, new PSF models could be make to match the lens and 
source star, and second round of fitting could be done with custom PSF models for the lens and 
source stars. 


An additional improvement in the method would be to fit more than two sources in the HST 
images in the vicinity of the target star. In the close-ups in the top-right and bottom two panels 
of Figure there is a faint star in the upper right corner. If this star were brighter or if the lens 
or source were much fainter, the PSF wings of this star could interfere with the lens and/or source 
star fits. The solution is then to also fit for that star. We expect to use this method for other 
targets that have been observed in Program GO-12541. 


Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the 
Space Telescope Science Institute (STScI), which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are 
associated with programs ^ I254I and I34I7. D.P.B.,A.B., and D.S. were supported by NASA 
through grants from the STScI and grant NASA-NNXI2AF54G. A.G. and B.S.G. were supported 
by NSF grant AST II0347 and by NASA grant NNXI2AB99G. S.D. is supported by the Strategic 
Priority Research Program- The Emergence of Gosmological Structures of the Ghinese Academy of 
Sciences (grant No. XDB09000000). The OGLE project has received funding from the European 
Research Gouncil under the European Gommunity’s Seventh Eramework Programme (PP7/2007- 
2013) / ERG grant agreement no. 246678 to AU. 
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